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Abstract:  This  paper  discusses  only  the  differences  between  the  original 
version  of  FASST  (Frankenstein  and  Koenig  2004a,  2004b)  and  the  new 
version.  This  report  is  intended  as  a  supplement  to  the  original  model 
documentation.  In  its  original  incarnation,  energy  and  mass  transport  as¬ 
sociated  with  water  vapor  in  the  soil  matrix  were  ignored.  The  author 
added  these  so  that  model  usage  could  be  expanded  to  include  biological 
investigations  yet  still  retain  its  initial  focus  of  soil  strength,  and  sensor 
performance  inputs.  Also  ignored  in  the  original  version  was  water  trans¬ 
port  due  to  soil  temperature  gradients.  In  determining  the  new  soil  tem¬ 
peratures  and  moistures,  the  original  model  first  achieved  convergence  in 
the  temperature  profile  followed  by  the  moisture  profile  at  a  given  time 
step.  The  new  version  of  FASST  solves  both  of  these  sets  of  equations  si¬ 
multaneously.  No  changes  have  been  made  to  the  equations  describing  the 
canopy  physical  state  except  to  allow  for  mixed  precipitation. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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Nomenclature 

/  subscript  indicating  low  foliage  (shrubs,  grass,  crops,  etc., 

vegetation)  terms 

g  subscript  indicating  ground  terms 

af,ag  shortwave  albedo  (o  -  l) 

avG  van  Genuchten  pressure  head  constant  (cm1) 
atj  numerical  calculation  constant 

fitj  numerical  calculation  constant 

Stj  numerical  calculation  constant 

At  time  step  (sec) 

Az  node  thickness  (m) 

sf ,  sg  longwave  emissivity  (o  -  l) 

sl  calculation  variable  (sl  =  sf  +sg~  £f£g ) 

(jkj  numerical  calculation  constant  (K) 

y  soil  water  surface  tension  (i g/s 2) 

y  =  75.6 -0.1425(7 -273. 15) -2.38x10^  (r -273. 15)2 

yd  bulk  dry  density  of  soil  ( kg/m3 ) 

yp  precipitation  density  ( kg/ms ) 

Yt  J  numerical  calculation  constant 

Y0  soil  water  surface  tension  at  25  °C  (71.89  g/s2) 

ij  mechanistic  enhancement  factor  for  vapor  flow 

Te  latent  heat  exchange  stability  correction  factor 

Th  sensible  heat  exchange  stability  correction  factor 
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r,„  momentum  stability  correction  factor 

k  thermal  conductivity  ( W/m-K) 

Ka  thermal  conductivity  of  dry  air  ( W/m-K) 

Ki  thermal  conductivity  of  ice  ( W/m-K) 

kv  thermal  conductivity  of  water  vapor  ( W/m-K) 

Kveg  vegetation  thermal  conductivity  ( Kveg  =  0.38  W/m-K) 
Kw  thermal  conductivity  of  water  ( W/m-K) 

vv  vertical  rate  of  water  vapor  flow  (m/s) 

vw  vertical  rate  of  water  flow  (m/s) 

W,i  numerical  calculation  constant 

Qi  volumetric  ice  content  (cm/cm^) 

Omax  maximum  soil  moisture  content  (m3/m3) 

6r  residual  soil  moisture  content  (m3/m3) 

Gv  volumetric  vapor  content  (cm3/cm3) 

Gw  volumetric  water  content  (cm3/cm3) 

pa  air  density  at  the  instrument  height  ( kg/m3 ) 

paf  density  of  air  in  the  atmosphere/foliage  interface 

Paf=Q-5{Pa+Pf)  (kg/m3) 

pag  density  of  air  at  the  foliage/ground  interface  ( kg/m 3) 

p  f  air  density  in  the  foliage  (kg/m3) 

Pi  density  of  ice  (kg/m3) 

Pv  density  of  water  vapor  (kg/m3) 

Pvs  saturated  water  vapor  density  (kg/m3) 

pw  density  of  water  (kg/m3) 

a  Stefan-Boltzman  constant  (5.6990-08  W/m2-K) 
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'F 

Cp 

Cp,a 

Cp.i 


CP,P 

Cp,s 

Cp.v 

Cp,veg 

Cp,w 

Ce 

c; 


Cf 

Ch 

cs 

'-/l 

c< 

'~hn 


cs 

hn 


Cn 


c,° 


foliage  fractional  coverage 

pressure  head  (m) 
specific  heat  of  the  soil  ( J/kg  -K) 
specific  heat  of  air  ( J/kg-K) 
specific  heat  of  ice  (-13.3  +  7.807^  J/kg-K) 
precipitation  specific  heat  ( J/kg-K) 
specific  heat  of  dry  soil  ( J/kg-K) 
specific  heat  of  water  vapor  ( J/kg-K) 
specific  heat  of  vegetation  (3500  <rr  J/kg-K) 

specific  heat  of  water  ( J/kg-K) 
empirical  coefficient  associated  with  Te 

bulk  transfer  coefficient  for  latent  heat  near  the  ground 


bulk  transfer  coefficient  near  the  ground  for  near-neutral 
conditions 

bulk  transfer  coefficient  for  turbulent  heat  in  the  foliage 
empirical  coefficient  associated  with  f;i 


bulk  transfer  coefficient  for  sensible  heat  near  the  ground 


bulk  transfer  coefficient  for  latent  heat  at  the  top  of  the  foliage 


for  near-neutral  conditions  C  fhn  = 


k!  In 


f  Z  -  Z,  A 

a  a 


V  J 


bulk  transfer  coefficient  for  sensible  heat  near  the  ground  for 
near-neutral  conditions 


empirical  coefficient  associated  with  fn 


bulk  transfer  coefficient  for  momentum  near  the  ground  for 
near-neutral  conditions 
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Cr  condensation  rate 

D  molecular  diffusivity  of  water  vapor  in  air  (m2/s) 

Dh  molecular  thermal  diffusion  coefficient  ( m2/s ) 

Dr,  Ds  vegetation  dripped  rain,  snow  (m/s) 

Dv  molecular  diffusion  coefficient  of  water  vapor  (m-/s) 

e0  windless  sensible  heat  correction  factor  (2.0  W/m2) 

E  evaporation  rate  (m/s) 

fi,f2,f3  variables  used  to  calculate  stomatal  resistance 
Ff  sum  of  energy  terms  at  the  atmosphere/foliage  interface  {W/m2) 

g  gravity  (9.81  m2/s) 

h  total  head  (m)  \  h  =  z-W  =  z-Pa/ pwg  ] 

hpond  head  due  to  water  collecting  on  the  surface  (m) 
hi, melt  head  due  to  melting  ice  (m) 
hs,melt  head  due  to  melting  snow  (m) 

Hf  sensible  heat  at  the  atmosphere/foliage  interface  {W/m2) 

H,  =  (e,  +l.lLAIpqfcraCrWqf)(Ta/  -  Tr) 

Hg  sensible  heat  at  the  foliage/ground  interface  {W/m2) 

Hs=Wp.tcpqlClWqf)(Tqr-Tg) 

Hr  relative  humidity  of  the  soil  (o  -  1) 

1 1  total  incoming  infrared  radiation  {W/m2) 

Ir,  Is  vegetation  intercepted  rain,  snow  (m/s) 

ij  total  incoming  solar  radiation 

k  von  Karmen’s  constant  {k  =  0.4) 

kth  soil  thermal  diffusivity  (kth  =  k/cp  m2/s) 

Km  pressure-driven  hydraulic  conductivity  (m/s) 

temperature-driven  hydraulic  conductivity  {m2/K-s) 


Kit 
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Ksat  saturated  hydraulic  conductivity  (m/s) 

Kvh  pressure-driven  vapor  conductivity  (m/s) 

Kvt  temperature-driven  vapor  conductivity  ( m2/K-s ) 

/  latent  heat  ( J/kg ) 

Ifus  latent  heat  of  fusion  (J/kg) 

hub  latent  heat  of  sublimation  (2.838eo6  J/kg) 

lw  latent  heat  of  vaporization  (J/kg) 

Lf  latent  heat  at  the  atmosphere/foliage  interface  ( W/m2) 

Lf  =  LAI pafCflWqfr  (qaf  —  qfsm) 

Lg  latent  heat  at  the  foliage/ground  interface  (W/m2) 

Ls=CselWafPag(qaf-qg) 

LAI  foliage  Leaf  Area  Index  (m2/m2) 

LAI  =  LAI^  +  (l  .0  -  0.00 1 6 [298.0  -  T,  ]2 )  [LAIm  -  LAImm  ] 

LAImin  minimum  foliage  Leaf  Area  Index  (rn2/m2) 

LAImax  maximum  foliage  Leaf  Area  Index  (m2/m2) 

mc  soil  clay  fraction  (o  -  l) 

mVG  van  Genuchten  constant  (mVG  =  l  -  l /tivg) 

M  molecular  weight  of  water  (0.018015  kg/mol) 

n  soil  porosity  (o  -  1) 

nVG  van  Genuchten  constant 

Pa  pressure  (Pa) 

Pf  precipitation  heat  at  the  atmosphere/foliage  interface  (W/m2) 

Pf  =-7P[  1.0-  exp  {-0.5  {LAI  +  &4/)}]  PcTp 

Pg  precipitation  heat  at  the  foliage/ground  interface  (W/m2) 

Pg  =  ~LP  [P,  ~  Interception  +  Drip}cp  Tp  (m/s) 

Pr  precipitation  rate  (m/s) 

qa  mixing  ratio  of  air  above  the  foliage 
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q,tf  mixing  ratio  of  the  air  at  the  foliage  interface 

c  _  [(l-crf)qa+CTf(0.3qa  +  0.6qf'Salr'  +  0.lqg'SatHR)] 

q“f  1  -  <rf  [0.6(1  -r")  + 0.1(1-^)] 

q/  mixing  ratio  of  air  at  the  top  of  the  foliage 

qg  mixing  ratio  of  air  at  the  ground  surface 

qfsat  saturated  foliage  mixing  ratio 

qg,sat  saturated  ground  mixing  ratio 

qtop  surface  water  flux  (m/s) 

ra  atmospheric  resistance  to  water  vapor  diffusion  ra  =  l/CfWaf 

(s/m) 

rce  turbulent  Prandtl  number  ( rce  =  0.71) 

rch  turbulent  Schmidt  number  (rch  =  0.63) 

r 

rs  stomatal  resistance  to  vapor  diffusion  rs  =  ff2f3  (s/m) 

rs,min  minimum  stomatal  resistance  to  vapor  diffusion  (s/m) 
r"  foliage  surface  wetness  factor  r"  =  r  (7;  +  rs )  1 

R  universal  gas  constant  for  air  (8.314  J/mol-K) 

Rib  bulk  Richardson  number 

SAI  Stem  area  index 

t  time  (s) 

T  temperature  (K) 

Ta  air  temperature  (K) 

Taf  air  temperature  in  the  foliage 

Taf  =  (l-*f)Ta+*f(03Ta+0.6Tf+0.1Tg)  (K) 

Tf,  Tg  foliage,  ground  surface  temperature  (K) 

Tp  precipitation  temperature  (K) 

77*  friction  velocity  (m/s) 
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Up  mass  precipitation  flux  ( kg/m2-s ) 

Vv  vapor  flow  rate  (m/s) 

Vw  water  flow  rate  (m/s) 

w  relative  soil  moisture 

W  wind  speed  at  the  instrument  height  Za  (m/s) 

Wqf  wind  speed  at  the  air/ foliage  interface 
Waf  =  0.83 <rfW'y[cl  +  (1  -  afW  (m/s) 

W'  adjusted  wind  speed  ( W'  =2.0  m/s  if  W is  below  2.0  m/s) 

z  depth  (m) 

zfo  foliage  roughness  length  (m) 

zg0  ground  roughness  length  (m) 

zho  sensible  heat  ground  roughness  length  (m) 

zqo  latent  heat  ground  roughness  length  (m) 

Za  shelter/instrument  height  for  air  temperature  (m) 

zd  zero  displacement  height  Zd  =  0.70 1Z°  975  (m) 

z f  vegetation  height  (m) 

zrh  height  above  the  ground  of  the  relative  humidity  measurement 

(m) 

Zu  height  of  the  measured  wind  speed  (m) 
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1.0  Introduction 

This  paper  discusses  only  the  differences  between  the  original  version  of 
FASST  (Frankenstein  and  Koenig  2004a,  2004b)  and  the  new  version. 

This  report  is  intended  as  a  supplement  to  the  original  model  documenta¬ 
tion.  In  its  original  incarnation,  energy  and  mass  transport  associated  with 
water  vapor  in  the  soil  matrix  were  ignored.  I  added  these  so  that  model 
usage  could  be  expanded  to  include  biological  investigations  yet  still  retain 
its  initial  focus  of  soil  strength,  and  sensor  performance  inputs.  Also  ig¬ 
nored  in  the  original  version  was  water  transport  due  to  soil  temperature 
gradients.  In  determining  the  new  soil  temperatures  and  moistures,  the 
original  model  first  achieved  convergence  in  the  temperature  profile  fol¬ 
lowed  by  the  moisture  profile  at  a  given  time  step.  The  new  version  of 
FASST  solves  both  of  these  sets  of  equations  simultaneously.  No  changes 
have  been  made  to  the  equations  describing  the  canopy  physical  state  ex¬ 
cept  to  allow  for  mixed  precipitation. 

To  begin,  section  2.0  provides  a  brief  synopsis  of  the  original  FASST  gov¬ 
erning  equations  and  solution  procedures  that  are  different  in  the  new 
version.  Readers  desiring  more  details  should  refer  to  Frankenstein  and 
Koenig  (2004a,  2004b).  Section  3.0  contains  the  new  version  of  the  equa¬ 
tions  and  solution  procedures.  A  complete  list  of  all  parameters  used  can 
be  found  in  the  beginning  of  this  report. 
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2.0  Original 

2.1  Temperature/Energy 

The  temperature  gradient  in  a  non-uniform  soil  layer  can  be  described  by 
the  one  dimensional  heat  flow  equation 

dT  lfus  p.  ddi  _  d  k  dT  f  ^  cp  w  QT 
dt  cp  pw  dt  dz[cp  dz  j  V"'  cp  dz 

where  T  is  the  temperature  (K),  t  is  time  (s),  k  is  the  thermal  conductivity 
(W/m-K),  cp,w  is  the  specific  heat  of  water  (J /kg  K),  cp  is  the  specific  heat 
of  the  soil  ( J/kg  K),  vw  is  the  vertical  rate  of  water  flow  (m/s),  lfus  is  the  la¬ 
tent  heat  of  fusion  (J/kg),  Oi  is  the  volumetric  ice  content  (cm^/cnP),  p\  is 
the  density  of  ice  (kg/ms),  pu,  is  the  density  of  water  (kg/ms),  and  z  is 
depth  (m)  measured  positive  downward  from  the  surface.  The  second  term 
on  the  left-hand  side  of  equation  (2.1)  represents  heat  lost/gained  due  to 
ice  formation/melting  and  the  terms  on  the  right-hand  side  incorporate 
temperature  changes  due  to  vertical  heat  conduction  and  water  flow,  re¬ 
spectively. 

If  low  vegetation  is  present,  the  atmosphere/foliage  energy  exchange  is 
given  as 


=  0  =  o-/  [/^  (!  -  af )  +  £fl\r  -  sfaT }  -  Pf 


(rf£f£g(J 


(TA-Tt)  +  Hf+Lf.  (2.2) 


2}  is  the  foliage  temperature  (K),  Tg  is  the  ground  temperature  (K),  eg 
is  the  ground  emissivity,  <7  is  Stefan-Boltzman  constant  (5.6990-08 
W/m2-K4),  if  and  if  are  the  total  incoming  solar  and  infrared  radiation 
(W/m2),  Hf,  Lf,  and  P/ are  the  sensible,  latent,  and  precipitation  heat  fluxes 
at  the  foliage  surface,  respectively  (W/m2),  and  =  sf  +  eg  -  sfsg .  The  foli¬ 
age  fractional  coverage  a/ ,  shortwave  albedo  a  f ,  and  emissivity  sf  are 
functions  of  the  vegetation  type  (high,  medium,  or  low)  and  season  (win¬ 
ter,  spring,  summer,  and  fall). 

The  energy  flux  exchange  at  the  foliage/ground  interface  as 
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F(Tg)  =  0  =  (l  -  af  )[^(1  -ag)  +  sjt  - e„oT4 


g  ir  g  g 


Gf£f£G  ,  . 

-Pg - -TAf) 


'l 


+  H  +L  +  k  — — +  / ,  —^-Az-v  c  T 

s  s  dz  fus  Pw  dt 


(2.3) 


where  (1  -  ay)  represents  the  radiant  fluxes  not  intercepted  by  the  vegeta¬ 
tion,  a„  is  the  shortwave  albedo  of  the  ground,  and  Hg,  Lg,  and  Pg  are  the 
sensible,  latent,  and  precipitation  heat  fluxes  at  the  ground  surface,  re¬ 
spectively  {W/m2),  and  A z  (m)  is  the  thickness  of  the  top  node.  Both 
and  are  a  function  of  the  soil  type  and  range  from  0.23  to  0.40  and 

0.92  to  0.97,  respectively.  The  third  term  in  the  second  row  of  equation 
(2.3)  takes  care  of  heat  conduction  to/from  the  surface  by  the  underlying 
ground,  depending  on  the  temperature  gradient.  This  is  followed  by  the 
heat  released/absorbed  by  the  soil  as  the  soil  moisture  melts/freezes.  Fi¬ 
nally,  the  last  term  represents  heat  that  is  advected  away  from/towards  the 
surface  as  a  result  of  the  vertical  movement  of  moisture. 


The  first  term  in  equations  (2.2)  and  (2.3)  represent  the  amount  of  solar, 
or  shortwave  radiation,  absorbed  by  the  surface.  The  second  term  is  the 
absorbed  incoming  longwave  radiation  while  the  third  term  is  the  emitted 
longwave  radiation.  The  sensible  and  latent  heat  fluxes  together  are  called 
the  turbulent  heat  fluxes  and  have  non-zero  values  in  the  presence  of 
wind.  The  precipitation  heat  represents  the  energy  needed  to  cool  or  heat 
any  snow  or  rain  that  falls  on  the  surface.  The  term  after  the  precipitation 
heat  flux  represents  the  radiative  exchange  between  the  foliage  and  ground 
surfaces.  In  equations  (2.2)  and  (2.3),  heat  that  is  transferred  to  the  sur¬ 
face  is  considered  positive. 

2.2  Moisture/Mass 

The  flow  of  water  ( vw )  through  a  porous  media  is  governed  by  Darcy’s  Law, 
which  states  that 


v  -K  —--K 

Vw  —  ^lh  ~  ^ Ih 

dz 


f  ,  FT  ^ 

-1  + - 

dz 


(2.4) 


where  Km  ( m/s )  is  the  pressure-driven  hydraulic  conductivity  and  h  (m) 
the  total  head  equals  the  elevation  head,  or  depth  (z),  minus  the  pressure 
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head  ( 'F ),  i.e.,  h  =  z-x¥  =  z-Pal pwg ,  Pa  (Pa)  is  pressure,  pw  ( kg/m3 )  is  the 
density  of  water,  and  g  (9.81  m2/s)  is  gravity.  For  unsaturated  soil,  'F  <  o. 

Also  governing  the  flow  of  moisture  through  a  soil  is  the  conservation  of 
mass,  which  states  that  the  time  rate  of  change  of  the  moisture  content  in 
a  given  volume  equals  the  net  gain/loss  of  fluid  in  the  volume,  i.e., 


39  dv 

w  _  _  w 

dt  dz 


Pw  & 


+  sources  -  losses 


(2.5) 


where  Qw  ( cms/cma )  is  the  volumetric  moisture  content.  Equation  (2.5)  as¬ 
sumes  that  changes  with  respect  to  time  in  the  soil  porosity  and  water  den¬ 
sity  are  negligible  compared  to  changes  in  the  soil  moisture  and  total  head. 
The  source  and  loss  terms  in  this  equation  account  for  occurrences  such  as 
runoff  and  plant  root  uptake. 

Equation  (2.5)  is  subject  to  the  following  flow  boundary  conditions  at  the 
surface  and  at  the  bottom  of  the  soil  column 

<ltop  =-E  +  C,  +  P+  (hpond  +  hi  melt  +  hs  melt )  /  A t  @  z  =  0 

Vbot  =  Kih  sin  (slope)  @z  =  zhot 


where  E  (m/s)  is  the  evaporation  rate,  Cr  (m/s)  is  the  condensation  rate, 

Pr  (m/s)  is  the  rate  of  precipitation,  hpond  (m)  is  the  head  due  to  water  col¬ 
lecting  on  the  surface,  and  h  i,melt  (m) and  h  s,melt  (m)  are  the  heads  due  to 
melting  ice  and  snow,  respectively,  and  At  (sec)  is  the  time  step.  If  the 
ground  is  sloped,  no  water  accumulates  and  any  water  that  falls  on  the 
surface,  but  does  not  infiltrate,  becomes  runoff. 

2.3  Hydraulic  Flow  Parameters 

The  pressure-driven  flow  parameters  are  unchanged  between  the  original 
and  new  versions  of  FASST.  Since  the  temperature-controlled  equivalent 
(see  section  3.3)  are  related,  I  include  them  here  for  reference  purposes. 
The  relationship  between  volumetric  moisture  content  and  pressure  head 
is  highly  nonlinear.  As  in  the  original  version  of  FASST,  following  the  work 
of  van  Genuchten  (1980) 
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0=6),+ 


6, . -6, 


(i+K^r)" 


(2.7) 


where  Or  is  the  residual  volumetric  water  content,  6max  is  the  maximum 
volumetric  water  content,  aVG  (cm-1)  is  a  constant  related  to  the  reciprocal 
of  the  bubbling  pressure  head,  nVG  is  a  constant  dependent  on  the  distribu¬ 
tion  of  pores,  and  mVG  =  1  -  1  /nVG. 


The  pressure-driven  hydraulic  conductivity,  Km,  is  modified  from  the 
original  version  of  FASST  to  take  into  account  the  decrease  in  flow  under 
freezing  conditions  and  is  given  as  (van  Genuchten  1980,  Niu  and  Yang 
2006) 


=K„(l+  W,  )V5 


w  =  ■ 


-  Or 

6 . -  6, 


(2.8) 


where  Ksat  (m/s)  is  the  saturated  hydraulic  conductivity. 

2.4  Turbulent  Energy  Drag  Coefficients 

The  bulk  transfer  coefficient  for  sensible  heat  C;f  (Deardorff  1978)  is  calcu¬ 
lated  using  the  bulk  transfer  coefficient  near  the  ground  Ct fH  (Koenig  1994) 
and  at  the  atmosphere/foliage  interface  C[n  (Balick  et  al.  1981)  for  near¬ 
neutral  stability  plus  a  sensible  heat  exchange  stability  correction  factor  rA 


cf  =  rt[(i-<7.)C'+c7.c/„] 


hn  1  ^  f^hn- 
2 


r^g  _ _ L 

^hn  r- 


ln(z.A‘) 


rch  + 


Hzjzj 
Hzjz f) 


(2.9) 


The  ground  roughness  length  zf  (m)  is  equal  to  0.001  m  for  all  soil  types 
and  0.0006  m  for  snow.  Since  Za  (m),  the  height  of  the  measured  air  tem¬ 
perature,  equals  Zu  (m),  the  height  of  the  measured  wind  speed,  Cfm  re¬ 
duces  to 
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cg  = 

hn 


(2.10) 


where  the  turbulent  Schmidt  number  rch  is  hardwired  in  the  code  as  0.63 
for  all  soil  types,  as  is  k,  von  Karmen’s  constant  (0.4).  The  term  Th  in  equa¬ 
tion  (2.9)  accounts  for  non-neutral  conditions  and  is  defined  as 


r 


h 


1-0  Rlh  <  0.0 

(l.0-16.0i?;7,)°5 

1.0  Rlb  =  0.0 

1.0 

1.0-5.07?/6  0.0  <  Rib  <  0.2 

R  JsZa{Taf-Tg) 

*  iTaf  +  Tg)wi 


(2.11) 


where  Rib  is  the  bulk  Richardson  number.  Similarly  for  the  latent  heat, 
r,  =  Th  and  Csen  follows  the  development  for  Cbn  with  Za  being  replaced 
with  Zrh  (m),  the  height  above  the  ground  of  the  relative  humidity  meas¬ 
urement,  rch  with  rce  the  turbulent  Prandtl  number  (0.71,  hardwired  in  the 
code). 

2.5  Precipitation  Flux 

In  the  original  version  of  FASST,  no  mixed  precipitation  was  allowed. 
Thus,  for  a  bare  surface  it  was  defined  as 


P  =U  c  T 

g  P  P>P  P 


U  =  —y  •  fallrate 


(2.12) 


with  Tp  the  precipitation  temperature  defined  as  the  air  wet-bulb  tempera¬ 
ture,  cp,p  as  either  the  specific  heat  of  water,  cp,w  (4217.7  J/kg-K)  or  ice,  cp,i 
( -13.3  +  7. 807;  J/kg-K)  depending  on  Tp,  Up  the  mass  precipitation  flux 
( kg/m2-s ),  and  yp  the  precipitation  density  ( kg/m.3 ).  Th e  fallrate  is  in 
(m/s). 
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2.6  Solution 

Since  the  solutions  to  the  temperature  and  moisture  balances  are  dis¬ 
cussed  in  detail  in  Frankenstein  and  Koenig  (2004a,  2004b),  I  present 
only  enough  material  so  that  comparisons  can  be  made  with  the  current 
procedures. 

2.6a  Temperature 

Equation  (2.1)  is  solved  using  a  modified  second-order  Crank-Nicholson 
approach.  Following  the  technique  presented  in  Hornbeck  (1975),  equa¬ 
tion  (2.1)  is  rewritten  as 


AS 

+ 

1 

AS 

hm  Pi 

1 

1 

1 _ 

1 

1 

+ 

1 _ 

At 

Cp  Pw 

At 

A z 

A z 

T  —  2T  +T  T  —2  T  +T 

±t+\,j+\  t+\J  ^  t  +  tyj+l  t,j  ^  t,j~\ 


(A*) 


(Az)2 


p,w 


w:t,j 


T  -T 

<j+i  i.j 

A z 


(2.13) 


where  kth  =  k/cp  and  the  subscripts  t  and  j  represent  time  and  depth,  re¬ 
spectively.  Combining  like  terms  and  rearranging  so  that  all  terms  involv¬ 
ing  Tj+ !  are  on  the  left-hand  side  of  the  equation,  equation  (2.13)  becomes 


-T 


Uj- 1 


+  YtjTt+\,j  +  Tt+\ j+i  ~ 

+  Tr,M,j  +  atJ  -  8tj)  +  Tt,j+ l(l  _  Pt.j  +  $tj)  +  Mtj  = 


where 


a,,j  =  “2 


'  (A^)2  ' 


kK,  (A0 


r,j=a,j- 2’  P,,j  =  2 


f  k  A 

lr 

V  J 


s'j=2\j 


f  .  A 
cp,w  Az 

\CP,,j  k<h,.j  J 


M,,j  =  atJ 


fus  Pi 
Cp  Pw 


(ML- 


(2.14) 


(2.15) 


For  more  detail  see  Frankenstein  and  Koenig  (2004a). 
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In  the  original  version  of  FASST,  the  soil  surface  and  foliage  temperatures 
were  solved  for  first.  In  order  to  solve  equations  (2.2)  and  (2.3),  we 
(Frankenstein  and  Koenig  2004a)  assumed  that 


t:= 


l  g,sat 


q 


g,sat 


(r/,+4(rX[VW, 
(fcLFilf-  [TATA 


tA7 iL 


(2.16) 


where  the  subscript  “t-i”  indicates  the  values  of  Tg  at  the  previous  time 
step.  7}4  and  qj,sat  are  also  represented  by  equation  (2.16),  substituting  the 
subscript  g  with/.  The  above  substitutions  allow  the  linearization  of  equa¬ 
tions  (2.2)  and  (2.3),  which  can  then  be  solved  simultaneously  for  7}  and 
Tg.  The  final  equations  are  thus  (Frankenstein  and  Koenig  2004b) 


c{  +  c{  Tg  +  c[  Tf  =  0; 
cf  +cfTg  +c\ \Tf  =  0. 


(2.17) 


2.6b  Moisture 

Equation  (2.5)  is  solved  numerically  using  an  explicit  scheme  such  that 


0  -0 

wi+i,j  w,j 

At 


( v  )  — (v  ) 

V  W/t+lJ+l  v  wj 


l+l  j 


A z. 


+  sources(j)  -  losses(  j) 


(2.18) 


where 


W,  =(*„),. 

(V«')y+1  —  1/2 

hj  =  Zj  -  KY  j  and  A zj  =  (z/M  -  zH )  /  2 .  The  change  in  soil  moisture  content 
due  to  changes  in  the  ice  content,  i.e.,  freezing/thawing,  is  incorporated 
into  the  source  and  sink  terms.  In  Equation  (2.5)  it  is  the  second  term 
on  the  right-hand  side.  Reexpressing  0w  in  terms  of  ip  using  equation  (2.7), 
equation  (2.18)  is  solved  for  ip  using  a  Newton-Raphson  technique. 


zj~zj-  1 


hj+i  ~  hj 
Z)+X~Zj 


(2.19) 
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3.0  New 


3.1  Temperature/Energy 

Equation  (2.1)  has  been  modified  to  include  temperature  changes  due  to 
vapor  flow.  It  is  now 


dcpT 

dt 


r  .  86 v 

fusPi  dt  +  WPW  dt 


8  f  8T  "N 

—  K- - VC  T-vc  (T - 273.15) -/ p  v 

V  dz  J 


dz 


+cpH,pw  ( sources  -  losses )  T 


(3-i) 


where  lw  is  the  latent  heat  of  vaporization  (J/kg),  6V  is  the  volumetric 
water  vapor  content  {cm3 / cm3),  cp,v  is  the  specific  heat  of  water  vapor 
( J/kg  K),  and  vv  is  the  vertical  rate  of  water  vapor  flow  (m/s).  Another 
change  from  the  original  formulation  is  that  z  (m)  is  now  measured  posi¬ 
tive  upwards  from  sea  level. 

The  calculation  of  the  specific  heat,  cp,  is  slightly  different  than  before  due 
to  the  incorporation  of  water  vapor  into  the  soil  matrix.  The  original  equa¬ 
tion  was 


=  {\-n)ydcP'S  +  6wp, 


Pp.w  +  (r 


-[0w+0i])PaCP,a+0iPiCp,i 


(3-2) 


It  is  now 

cp=(\-  n)ydcps  +  6wpwcpw  +  (n  -  [6W  +  6>  +  6^)pacpa  +  6ipicpi  +  6vprcpv .  (3.3) 

See  the  list  of  variables  at  the  beginning  of  the  report  for  a  description  of 
all  of  the  parameters.  While  the  calculation  of  the  thermal  conductivity,  k, 
is  the  same  as  in  the  original  version  of  FASST,  the  calculation  of  the  indi¬ 
vidual  components  is  different.  Comparison  of  the  original  and  new  values 
of  the  thermal  conductivity,  specific  heat,  and  density  for  water,  air,  ice, 
and  water  vapor  are  listed  in  Table  1.  Unlike  in  the  original  version  of 
FASST  where  we  assumed  that  the  specific  heat  and  thermal  conductivity 
of  vegetation  were  negligible,  I  now  assume  that  they  are  cp,veg  =  3500 
J/kg  K and  Kveg  =  0.38  W/m-K  (Moore  and  Fisch  1986). 
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Table  l.  Physical  properties  of  water,  ice,  air,  and  water  vapor. 


Variable 

Old 

New 

Cp,w 

T(°C) 

(J/kg-K) 

4217.7 

-10-6775  +  10-4T4  -  6  X10-3T3  +  lO-'T2  -  2.9 7  +  4216.9 

Hillel  (1998)  4216.9  <  cP;W  <  31500.0 

Cp,a 

T(°C) 

(J/kg-K) 

1005.6 

1.9327  x  10-9T4  -  7.9999  x  10-773  +  1.1407  x  10-371= 

-0.4489T  +  1057-5 

Ierardi  (1999)  1000  <  cp,a  <  1250 

(J/kg-K) 

-13.3  +  7.807; 

-13.3  +  7-8o7a 

Jordan  (1991)  1389  <  cp,i  <  2050 

Cp,v 

(J/kg-K) 

2.3888  x  io-874  -  6.5129  x  io-573  +  6.6178  x  10 ~2T2 
-2.9086  x  1047+  6625.5 

Ierardi  (1999)  2000  <  cp,v  <  3260 

pw 

(kg/m?) 

1000.0 

-3.0  x  10-774  +  7.0  x  10-573  -  9.92  x  10-373 
+8.666  x  io~27+  999.81 

Hillel  (1998)  962  <  p„  <  1000 

Pa 

(kg/mp 

3.48x10  3 

36o.778i97-10°336 

Ierardi  (1999)  0.948  <  pa  <  2.05 

Pi 

T(°C) 

(kg/mp 

916.5 

-2.0  x  io-1076  -  7.0  x  io-873  -  1.0  x  10-574 
-7.0  x  10-473  -  2.37  x  io-273  -  4.36  x  10-473  +  916.2 
Engineering  ToolBox  (2005)  916.2  <  pi  <  925.7 

Pv 

(kg/mp 

4.192  x  io-i274  -  1.25128  x  io-873  +  1.45079  x  10-572 
-8.12253  x  10-37+  2.17634 

Ierardi  (1999)  0.60  <  pv  <  1.14 

Kw 

(W/m-K) 

1.8  xlCT3  T  -0.0787 

1.8  x  10-37-  0.0787 

Farouki  (1981)  0.58  <  kw  <  0.75 

Ka 

(W/m-K) 

1.5207  x  10-1173  -  4.8574  x  io-873  +  1.0184  x  10-47 

-3-9333  x  10-4 

Ierardi  (1999)  0.0159  <  ica  <  0.1000 

Ki 

(W/m-K) 

488. 19I-1  +0.4685 

4.0  x  10-773  +  1.0  x  10-472  -  6.9  x  io-37+  2.2174 

Engineering  ToolBox  (2005)  2.2174  ^  k;  <  3.4800 

Kv 

(W/m-K) 

8.3154  x  10-57-  7.4556  x  10-3 

Ierardi  (1999)  0.00695  ^  kv  <  0.02360 

The  new  surface  boundary  conditions  at  the  soil  and  vegetation  surfaces 


are 
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Ff=  0  =  <jf 


li<X-af)  +  efli-efoT}-Pf 


Gf£f£G  . 

+  f  f  g  (T*-Tf) 


+H f+Lf+  G fKveg 


dTf 

dz 


(34) 


F(  T„  )  =  0  =  ( 1  -  <7, )  [//  (1  -  a, )  +  sJi  -  st  < 


Gf£f£0G  ,  , 

-P„ - L-Ljl-(T,t-T,f) 


+H,+L,+l,..^-^LAz-lwpw^  +  K^- ~(vwcr  +  y,c„ )(rs  - 273 . 1 5) 


lg  '  '  - Jus 


pw  dt  dt  dz 

-lwPwvv  +  cp,wPw  ( sources  -  losses)T  . 


(3-5) 


Unlike  in  the  original  version  of  FASST  where  the  bottom  temperature  was 
allowed  to  float,  I  assume  a  constant  deep-earth  heat  flux  for  the  bottom 
boundary  condition  of  75  mW/m2. 

3.2  Moisture/Mass 

The  water  flow  rate  ( vw )  is  now 


-K  —-KT  —--K 
^ Ih  ~  ^IT  ^  ^lh 

oz  oz 


A  dy 
1  +  — 
v  dz 


-K 


IT 


dT_ 

dz 


(3-6) 


where  Kit  is  the  temperature-dependent  hydraulic  gradient  ( m2/K-s ).  The 
new  moisture  governing  equation  takes  into  account  water  flow  due  to 
temperature  gradients,  unlike  the  old  version,  which  considered  only  flow 
due  to  gravity  and  pressure  gradients.  As  with  the  temperature  governing 
equations,  vapor  fluxes  are  also  now  accounted  for.  The  new  mass  balance 
equations  is 

d#w  +  d0v  |  a  d6 \  _  d 

dt  dt  pw  dt  dz 

(3-7) 


(AT  dT 

K„1-^(K„,+K„)  +  —(Kit+K,t) 


+  sources  -  losses 


where  Kvr  is  the  temperature-dependent  vapor  gradient  ( m2/Ks ),  and  Kvh 
is  the  pressure-dependent  vapor  conductivity  (m/s).  The  surface  moisture 
boundary  condition  remains  unchanged  from  the  original  version  of 
FASST. 


ERDC/CRREL  TR-08-7 


12 


Differences  between  equations  (2.4)  and  (3.6)  can  be  explained  by  the 
change  in  the  vertical  reference  (z  positive  downward  from  the  surface 
previously  versus  z  positive  upward  from  sea  level  currently)  and  the  in¬ 
clusion  of  temperature  gradient  flow. 

3.3  Hydraulic  Flow  Parameters 

The  temperature-dependent  hydraulic  gradient,  Kit,  is  given  as  (Hansson 
et  al.  2004,  Noborio  et  al.  1996) 


K 


IT 


=  K„ 


f 

7' ¥ 

v 


1  dy^ 

Yo  dT  , 


(3-8) 


y  =  75 .6  -  0. 1 425  ( T  -  273 . 1 5)  -  2 .3 8x1  (T4  ( T  -  273 . 1 5)2  is  the  soil  water  sur¬ 
face  tension  ( g/s 2),  y0  =  y  (25  °Q  =  71.89  g/s2,  and  Km,  the  pressure-driven 
hydraulic  conductivity,  is  defined  in  equation  (2.8). 


The  water  vapor  flow  rate  ( vv )  is  given  by  (Fayer,  2000) 


V  —K  d-^-K  <H- 

Vv  v/i  -  YLvT 

oz  oz 


where  (Fayer  2000,  Noborio  et  al.  1996,  Nassar  and  Horton  1989) 

D  =  0.66[n -(0„  +  6t)] 


_  DpvsMg 

-<v Vh  -  r™  nR 


2.12x10^ 


'  t  5 

2 

L  273.15 J 

Kvt  = 


pwRT 

DrjHR  dpvs 
A,  dT 


7=9.5  +  66 w  - 8.5  exp  j -[(l  +  2.6m;112 )  Gw 

'FMg 


Hr  =  exp 


RT 


(3-9) 


(3-10) 


D  is  the  molecular  diffusivity  of  water  vapor  in  air  ( m2/s ),  n  is  the  porosity, 
77  is  the  mechanistic  enhancement  factor,  mc  is  the  clay  fraction,  Hr  is  the 
soil  relative  humidity,  M  is  the  molecular  weight  of  water  (0.018015 
kg/mol),  R  is  the  universal  gas  constant  (8.314  J/mol-K),  and  pvs  is  the 
saturated  water  vapor  density  (kg/m3). 
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3.4  Turbulent  Energy  Drag  Coefficients 

An  effort  was  made  to  move  away  from  a  more  empirical  formulation  of 
the  turbulent  energy  terms  and  instead  to  base  them  on  the  more  widely 
accepted  Monin-Obukhov  similarity  theory.  Unfortunately,  this  requires 
iterating  for  a  solution,  which  is  numerically  cost-prohibitive.  In  order  to 
avoid  this,  I  adopted  the  method  of  Mascart  et  al.  (1995)  and  Louis  (1979) 
The  bulk  transfer  coefficients  Cfm  and  C8en  are  now  defined  as 


rg  = _ 

^ hn 

In 

(3-n) 

<2=- 
in  | 

where  zho  (m)  is  the  sensible  heat  roughness  length  and  zqo  (m)  is  the  mois¬ 
ture  roughness  length  and  all  other  terms  are  described  in  Section  2.3.  The 
roughness  lengths  are  calculated  using  (Jacobson  2005) 


h  . 

z„  =  S 


max 


zq  =  < 


max 


(  D  x 

0-9  zl,-± 

y  KU*  j 


f  D  x 

0-9zf,-^ 

y  Ku*  j 


ut  =  0.0 
iu  ^  0.0 

ih  =  0.0 
Ut  ^  0.0 


A 


K . 


PaCr, 


a  p,w 

if  Dh  <  0.0,  Dh  =  Dv 


D, 


2.1  lxl  O'5 


(  T  ^ 

a 

1.94 

1013.25^ 

V 273.15  j 

p 

\  ra  J 

(3-12) 


with  (Louis  1979) 
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u;  =  c°  w2 r 

*  ng  m 


C°  = 

ng 


In 


f  7  /  \ 


V  /  z0  J 


r_  = 


l- 


9.4  R, 


1  +  Cm  Rib 


1 

I  +4.7/?,:, 


1 1  /  2 


*  0 


^  >  0 


CL  =  7.4C”  9.4  1^- 

\/  g 

Z0 


(3-13) 


where  f?a,  the  bulk  Richardson  number,  is  defined  as 


R 


ib 


0.21 

.  h/+7;K 


»V<o.oi 


»V*o.oi 


and  ra/(K)  and  Wa/  (m/s)  are  the  temperature  and  wind  speed  at  the  foli¬ 
age/ground  interface  respectively  (Deardorff  1978).  The  sensible  heat  ex¬ 
change  stability  correction  factor  term  (170  in  equation  (2.9)  that  accounts 
for  non-neutral  conditions  is  defined  as  (Mascart  et  al.  1995) 


r/;  =1 
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V  /  z0  J 


v  /  zoy 
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1 

1  + 4.7ft, 


|l/2 


R,b  *  0 


Rib  >  0 


(3-14) 


The  latent  heat  exchange  stability  correction  facter  term,  Te,  is  the  same  as 
Th,  replacing  zhQ  with  zqa  and  Ch  with  Ce.  Ch  is  defined  as  (Mascart  et  al. 
1995) 
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C  =9  4CC 


Vz°  J 


Cl  =  3.2 165  +  4.343 1//  +  0.536//2  - 0.078 1//3 
ph  =  0.5802-0.1571//  +  0.0327/r-0.0026//3 


//  =  In 


(3-15) 


Ce  is  similarly  calculated  by  substituting  z*  with  zqo . 

3.5  Precipitation  Flux 

In  the  new  version  of  FASST,  mixed  precipitation  is  allowed.  Thus,  for  a 
bare  surface  it  is  defined  as 

Pg=~ (PWCP,W  '  [rain  fallrate  - Ir  +  D, .]  +  ypcpi  ■  [snow  fallrate -  Is  +  Ds ]) Tp 


(3-16) 


with  Tp  the  precipitation  temperature  equal  to  the  air  wet-bulb  tempera¬ 
ture,  Ir,  Is  are  respectively  the  vegetation  intercepted  rain  and  snow  (m/s), 
and  Dr,  Ds  are  the  vegetation  dripped  rain  and  snow  (m/s),  respectively. 

3.6  Solution 

Unlike  the  original  version  of  FASST  where  the  temperature  profile  was 
solved  for  before  the  moisture  at  a  given  time  step,  in  the  new  version  both 
profiles  are  iterated  for  simultaneously.  Another  change  from  the  original 
version  is  that,  since  the  orientation  of  the  z-axis  is  reversed,  node  1  is  now 
at  the  bottom  of  the  soil  instead  of  the  top.  This  allows  the  expansion  of 
the  solution  matrix  if  snow  is  present. 

3.6a  Temperature 

Following  the  technique  outlined  in  Celia  et  al.  (1990)  and  Hansson  et  al. 
(2004)  the  individual  terms  in  the  left-hand  side  of  equation  (3.1)  are  re¬ 
written  as 
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where  8]  =  r“  -  7  ,  fl  and  the  subscript  is  depth  and  the  superscript 
“k”  is  the  iteration  level.  The  right-hand  side  of  equation  (3.1)  for  the  inte¬ 
rior  nodes  becomes 
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Expanding  only  the  q/+ 1/2,2+ 2/2  term  since  the  qj-i/a,t+i/2  term  is  similar 
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If  vegetation  is  present,  then  equation  (3.20)  for  the  vegetation  (j  =  nnodes 
+  1)  and  surface  (j  =  nnodes)  nodes  is  written  respectively  as 
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(3-24) 


Substituting  equations  (3-i9)-(3.24)  into  equation  (3.1),  assuming  (2.16) 
is  still  valid  and  rearranging  so  that  all  8T  terms  are  on  the  left-hand  side, 
the  final  temperature  equation  is 


A)  8j_x  +  B.  Sj  +  Cj  Sj+1  =  D.  . 


(3-25) 


A  Newton- Raphson  technique  is  used  to  solve  the  system  of  linear  equa¬ 
tions  as  with  the  original  version  of  FASST. 

3.6b  Moisture 

Similar  to  the  solution  technique  for  the  temperature  equation,  I  follow  the 
technique  outlined  in  Celia  et  al.  (1990).  Equation  (3.7)  is  thus 
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where 


S0_ 

dt 


80 w  80v  pt  80j  , 

■  +  — L  +  — — — ,  v  =  v  +  v  and  v 


dt  dt  pw  dt 


j±l/2,t+l/2  n  y  j±l/2,t+l  Vj±±/2. 


At  the  surface,  Vj+i/2,t  =  qtop  while  at  the  bottom,  =  qbot,  both  of  which 
are  given  in  equation  (2.6).  Expanding  the  left-hand  side  of  equation 
(3-26) 
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(3-27) 


where  “k”  is  the  iteration  level  and  A'1'  =  t+1 .  The  [vw]y+1/2  r+1/2  term 

in  the  right-hand  side  of  equation  (3.26)  becomes 
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Except  for  the  “+  1”  terms  in  equation  (3.28),  the  expression  for 
[vv].+1/2 M| /2 is  the  same.  The  same  principle  is  applied  to  the  [v]  /+1/2 

term  in  equation  (3.26).  Substituting  equations  (3.27)  and  (3.28)  into 
equation  (3.26)  and  rearranging  so  that  all  AT  terms  are  on  the  left-hand 
side,  the  final  moisture  equation  is 
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AJSJ_1+BJSJ  +  CJSJ+1=DJ 


3-29 


As  with  the  original  version  of  FASST,  a  Newton- Raphson  technique  is 
used  to  solve  the  system  of  linear  equations. 
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